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Abstract 

We consider finite sequences s € D n where D is a commutative, unital, integral 
^ \ domain. We prove three sets of identities (possibly with repetitions), each involv- 

ing 2n polynomials associated to s. The right-hand side of these identities is a 
recursively-defined (non-zero) 'product-of-discrepancies'. There are implied itera- 
tive algorithms (of quadratic complexity) for the left-hand side coefficients; when 
qs ' the ground domain is factorial, the identities are in effect Bezout identities, 

r**"- . We give a number of applications. Firstly, a new (quadratic) algorithm to 

compute Bezout coefficients over a field, which compares favourably with the ex- 
tended Euclidean algorithm. We show that the successive output polynomials of 



(N 
> 



■ the Berlekamp-Massey algorithm either coincide or are relatively prime. A second 

application concerns sequences with perfect linear complexity profile. We give new 
characterisations of them in terms of minimal polynomials and a simpler proof of a 
theorem of the characterisation of binary sequences with perfect linear complexity 
^ . due to Wang and Massey using the third set of identities. 

Another application is to annihilating polynomials which do not vanish at zero 
and have minimal degree. We simplify and extend an algorithm of Salagean to 
sequences over D. First we prove a lower bound lemma which was stated without 
proof for sequences over a field. This lemma and an easy extension of the author's 
minimal polynomial algorithm yields such annihilators. The first set of identities 
allows us to remove a test. In fact, we compute minimal realisations and give 
corresponding identities. We also construct these annihilators by extending the 
sequence by one term, give corresponding minimal polynomial identities and prove 
a characterisation (stated without proof by Salagean over a field). In the Appendix, 
we give an alternative proof of the lower bound lemma using reciprocal annihilators 
and apply it to the complexity of reverse sequences. This gives a new proof of a 
theorem of Imamura and Yoshida on the linear complexity of reverse sequences, 
initially proved using Hankel matrices over a field and now valid for sequences over 
a factorial domain. 
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1 Introduction 

1.1 Overview 

Let F be a field and (u, 112) G ¥[x} 2 not both zero. It is well-known that there are Bezout 
coefficients (/, f?) G ¥[x] 2 such that f-u+f2-U2 = gcd(w, 112), known as a Bezout£j] identity. 
(Of course, it easy to check that if (/, ^2) are Bezout coefficients for (m,m 2 ), then so are 
(/) f-i) +h ■ (—U2, u)/ gcd(w, u 2 ) for any h G ¥[x}.) As far as we know, Bezout coefficients 
are always computed using the extended Euclidean algorithm. 

In this paper, we are interested in the case when (u, U2) arises from a finite sequence. 
For example, this pair could be the pair of polynomials used to find a minimal one, the 
successive minimal polynomials or the rational function approximating the 'generating 
function' of the sequence. All of these are typically relatively prime polynomials. These 
pairs are all computed using the minimal polynomial algorithm due to the author [T2] . 
|14j ; for a succinct overview, see [IT] . A natural question is whether there are associ- 
ated minimal polynomial identities. We answer this in the affirmative for each pair of 
polynomials and give some applications. 

Thus let s = (si, . . . , s n ) G D n denote a sequence of n terms from a commutative unital 
integral domain D, a 'domain' for short. We prove theorems giving n minimal polynomial 
identities (possibly with repetitions) associated to s, Theorems 16.41 and 16.81 It turns out 
the right-hand sides are equal to the same recursively-defined 'product-of-discrepancies' 
V n G D \ {0} of Definition 16.31 There is an implied iterative algorithm of quadratic 
complexity, Algorithm 16.131 (When the ground domain is factorial, the identities are in 
effect Bezout identities as the minimal polynomials are relatively prime.) 

We give a number of applications: the successive 'connection polynomials' of the 
Berlekamp-Massey algorithm (as developed in [10]) either coincide or are relatively prime 
when D is factorial, Corollary 16.21 and a new (quadratic) algorithm to compute Bezout 
coefficients, which compares favourably with the extended Euclidean algorithm; see Sec- 
tion 7. Another application concerns sequences with perfect linear complexity profile. In 
Section 8 we give new characterisations of them in terms of minimal polynomials and 
a simpler proof of a theorem of Wang and Massey characterising binary sequences with 
perfect linear complexity, [20J. 

Section 9 gives applications to annihilating polynomials which do not vanish at zero. 
We simplify and extend to D results of [18], which used definitions of a linear recurring 
sequence, annihilating and minimal polynomial for finite sequences over F (with first 
term sq) equivalent to the non-standard definitions of p2], [14^1 . First we establish a 
lower bound lemma which was stated without proof in [15] . The resulting corollary and 

1 Etienne Bezout (1730-1783); according to [HI p. 115], the integer case is due to Claude Bachet de 
Meziriac (1581-1638): "Problemes plaisans et delectables qui se font par les nombres", 1624. 
2 We note that [12] and [14] were cited in [15] . 
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algorithm do not require any characterisations of minimal polynomials. In particular, the 
derivation of Algorithm 19.61 is similar to [6] and [121 Algorithm 4.2]. Furthermore, our 
first identity allows us to eliminate a test used in [18, Algorithm 3.2]. We also construct 
annihilators that do not vanish at zero by extending the given sequence by one term and 
give corresponding minimal polynomial identities, Theorem 19.91 We characterise these 
annihilators of minimal degree in Theorem 19.13} a version over F was stated without 
proof in [IB] . 

In the Appendix (Section 10), we give an alternative proof of Lemma [9.31 using recip- 
rocal annihilators and apply it to the complexity of reverse sequences. This gives a new 
proof of [4, Theorem 3, p. 149] on the linear complexity of reverse sequences, initially 
proved using Hankel matrices over a field and now valid when D is factorial. 

A preliminary version of this work was presented in May 2010 at Equipe SECRET, 
Centre de Recherche, INRIA Paris- Rocquencourt, whom the author thanks for their hos- 
pitality. 

1.2 In More Detail 

Our approach is based on the recursive minimal polynomial theorem of [H]] which is 
relative to any e E D. This gives the next minimal polynomial /i^ using a previous 
one /i^" --1 ) and a //( n_1 ); informally, /i^™ -1 ) is a 'prejump' polynomial of (si, . . . , s n _i), 
i.e. a minimal polynomial immediately preceding a jump in the linear complexity L„_i = 
deg(/i ( - n_1 ^) — we are ignoring the initial cases. We prove inductively that /i*-" - -' and /i'*-™* 1 
are always relatively prime (when D is factorial) and exhibit polynomials and /g 
such that for any domain D 

fin) . ^n) + fin) . ^n) = ^ 

We apply this to obtain an identity for /i*™ -1 ) and Proposition 16.61 

Before discussing our third identity, we need to define another polynomial /4"^ derived 
from // n ) and s; /i^ G D[x] is the polynomial part of the Laurent product 

■ (s n x~ n H h SiaT 1 ). 

In an earlier paper [12] , we inductively constructed jffl = (fi^\ fi^) £ -D[s] 2 for 1 < i < n 
(possibly with repetitions) such that 

(i) deg(/4f ) ) < deg(/i«) 

(ii) fi« ■ (six-* + ■■■ + siar 1 ) - nf = mod ( x de g ( M «)-i-i) 

(iii) [x^ is a solution to (ii) of minimal degree. 

See Algorithm 16.131 which does not use scalar or polynomial division. We called the pair 
minimal realisation' of s. If /i ( - n - ) is monic, this gives a rational approximation of 
Six^ 1 + • • • + s n x~ n with denominator of minimal degree. When D is factorial, minimality 
easily implies that fi^ and /4 are relatively prime, so it is natural to ask what their 
'coefficients' are. Our third identity is 
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Algorithm 16.131 already computes these coefficients, so they can be read off at the end 
of the n th iteration at no extra cost. The degrees of the coefficients are bounded above 
by L n _i. We also show that if D is a field, these coeffients are valid for any minimal 
realisation of s, Corollary 16 .101 Secondly, our coefficients agree with the coefficients found 
by the extended Euclidean algorithm (up to a scalar), Proposition 16.111 

We conclude this Introduction by outlining our algorithm for obtaining the Bezout 
identity for (w, 1*2) £ D[x} 2 where < deg(w 2 ) < deg(w) = d, u is monic and D is a 
principal ideal domain. The steps are: 

(i) 2d subtractions in D[x~ l ,x] to obtain s = U2/u(mod x~ 2d ^ 1 ) 

(ii) Algorithm 16.131 for the Bezout coefficients associated to s 

(iii) two multiplications in D[x) to obtain V2d • gcd(w, M2). 

This means that we can use Algorithm 16.131 whenever the extended Euclidean al- 
gorithm is applied e.g. to compute Bezout coefficients in F[x,y] for example. Bezout 
identities also arise in the context of the key equation of Coding Theory, see e.g. [7], 
|13j . [T5] and it seems likely that the algorithm of [21 Section 4] is closely related to the 
minimal realisation algorithm of [T2]; see also Section 8]. 

2 The Recursive Minimal Polynomial Theorem 
2.1 Notation 

Our notation is continued from [10J : N = {1, 2, . . .}, n G N and D is a commutative, unital 
integral domain with 1 ^ 0, or 'domain' for short. For any set S containing 0, S x = S\{0}. 
We say that / G _D[x] x is monic if its leading term is 1 and the constant term of / is /o- 
The reciprocal of is and for / G D[x] x , its reciprocal is f*(x) = x deg ^f(x' 1 ). We 
often write / = x e g + h for f(x) = x e g(x) + h(x), where eGff and g,h G D[x\. 

For / G D[x], f denotes an element of Z)[x] 2 with first component / and a certain 
component fi (depending on a sequence) to be specified. We regard D[ clS cL D\x\- 
module by multiplication on each component; e.g. for e G N, A G D and /, /' G D[x] 

x e J + A-j' = (x e f + A • f',x e f 2 + A • f 2 ). 

We will also extend deg to D[x} 2 via deg(/) = (deg(/), deg(/ 2 )) and write gcd/ for 
gcd(/,/ 2 ). If7=(/,/ 2 ),welet 

/=(-/a,/) 

We also use • for the inner product f -u — (fu, Z^"^)- 

We write : Z)[[x _1 ,x] for the domain of Laurent series in x~ l over D and let v : 
/^[[x" 1 ,^] — > Z be the exponential valuation. We also use v denote its restriction to 
D[x~ l ,x]. Thus v : — > Z is v(g) = max{i : gi ^ 0}, which coincides with the 
degree function on D[x}. We have v(gg') = v(g) + v(g'), v(g + g') < m&x{v(g), v(g')} and 
y(9 + 9 1 ) = max{v(^), v(g')} if v(g) ^ v(g)'. 
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2.2 Sequences 

Our non-standard approach to finite sequences was based on regarding linear recurring 
sequences s : N — > D as the torsion submodule of a natural .D[x]-module, [121 Section 2], 
[TU Section 3]. We begin with Laurent series in x , -D[[x -1 ] as standard D[[x _1 ]-module 
i.e. acting on itself via multiplication. This also makes Z)[[x _1 ] into a D[x]-module. Let 
s = J2j>i s j x ~ J '• Then for / G D[x], put 

i>i 

One checks that o makes x _1 -D[[x -1 ]] into a -D[x]-module. 

Definition 2.1 (Cf. fl?^ f^, JT8$) We say that s satisfies a linear recurrence relation 
if it is a torsion element i.e. if Ann(s) = {/ G D[x] : / o s = 0} {0} or if for some 
f G D[x] x (/ • s)<z-j = for d — j < —I where d = deg(/) i.e. 

foSj-d H h /dSj = /or d + 1 < j. 

When fd = 1, we can write Sj = —{foSj-d + • — h fd-iSj-i) for j > d + 1 and s is a linear 
recurring sequence. 

Now for n G N and s = (si, . . . , s n ) G -D n , let s G -D[x] be = six + • ■ ■ + s n a; n . We 
will also abbreviate = six' 1 + • • ■ + s n x~ n to s, so that Sj = s _j for — n < j < —1. 

In the following definition, multiplication of / G -D[x] and s G -D[x _1 ] is in the domain of 
.D-Laurent polynomials 

Definition 2.2 (Annihilator, annihilating polynomial) ( JJE. Definition 2.7, Propo- 
sition 2.8], Cf. fl^, fTBJ)) If s G D n , then f G D[x] is an annihilator (or a characteristic 
polynomial) of s if f = or (/ • = /or — n < d — j < —1 i.e. 

foSj-d H h /rfSj = for d + 1 < j < n 

where d = deg(/) > ; written f G Ann(s). 

Any polynomial of degree at least n annihilates s, vacuously. For 1 < % < n, we write 
S W for (si, . . . , Si). The following definition is a functional version of [12, Definition 2.10]. 

Definition 2.3 (Discrepancy Function) We define A : D[x] x x D n D by 

A (/> s ) = (/-s)deg(/)-n- 

Thus A(/, s) = J2k=o fk s n-d+k where d = deg(/). Clearly for n > 2, / G Ann(s) x if and 
only if / G Ann(s < - n_1 - ) ) >< and A(s,f) = 0. If s G -D n is understood, we write A n (/) for 
A(/, s); if / is also understood, we simply write A n . It is elementary that if 1 < i < n — 1, 
then (/ • sl!l)de g (/)-i = (/ • s)de g (/)-i- 
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3 Minimal Polynomials 

Definition 3.1 (Minimal Polynomial) f/HJ Definition 3.1], Cf. /]]/, Hfj We say 
that f G Ann(s) is a minimal polynomial of s G D n if 

deg(/) = min{deg(o) : a G Ann(s) x } 

and let MP(s) denote the set of minimal polynomials of s. 

As any / G of degree at least n annihilates s G D n , MP(s) 7^ 0. We do not require 
minimal polynomials to be monic. For any d G D x , d G MP(0, . . . , 0); if s\ ^ and 
deg(/) = 1 then / G MP((si)) since D has no zero divisors. 
The linear complexity function L : D n — > {0} U N is 

L(s) = deg(/) where / G MP(s). 

We will also write L n for L(s) when s is understood and similarly Lj = L(s^) for 1 < 
j < n. For fixed s, L is clearly a non-decreasing function of i. 

3.1 Exponents 

Definition 3.2 (Exponent Function) Forn > 1, let the exponent function e n : D[x] x x 
N — ?• Z &e given fry 

e n (/)=n + l-2deg(/). 

The following lemma is the annihilator analogue of (6J Lemma 1]. 

Lemma 3.3 (JZB £emma 5.K/j Letn>2, f e Ann(s (n " 1) ) x and A n (/) ^ 0. 
(%) For any g G Ann(s) x ; deg(o) > n - deg(/) = e n _i(/) + deg(/). 
If h E MP(s) tfien deg(a) > max{e n _i(/i), 0} + deg(/t). 

4 A Recursive Minimal Polynomial Function 

We will define a recursive minimal polynomial function /x : -D n — > D[x\. But first we need 
the following function (which assumes that fi : D n ~ 1 — > D[x] has been defined). We also 
set A = 1. 

4.1 The Index Function 

Definition 4.1 (Index Function) Let n > 1 and s G D n . We set fi^ = 1 (so that 
Ai = A^// )) = sj and e = 1. Suppose that for 1 < j < n - I, G MP(s^O) has 
discrepancy Aj+i and exponent ej. We define the index function 

':{0,...,n}^{-l,n-l} 

by 0' = — 1 and /or 1 < j < n — 1 

•/ _ f U ~ */ A i = or (Aj ^ and e,-_i < 0) 
~~ \ j - 1 if Aj and e 3 -_i > 0. 
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4.2 The Recursive Theorem 

We recall the recursive function ft : D n — >■ D[x] x of [ID] : for all s G -D n , /i(s) G MP(s). 
When s is understood, we will write for fj,(s^). 

Definition 4.2 (Basis of the Recursion) Recall that 0' = — 1 and A = 1. Let e G D 

be arbitrary but fixed and s G D n . We put n^~ l > = /i(s, —1) = e and ^S Q > = jj(s, 0) = 1. 

Thus the exponent of /i^ is eo = 1 and Ai = sj. We also recall some notation from [TO] : 

(i) //' = fi o ' (where o denotes composition) 

(ii) 17 = deg o// 

(iii) A' = A o (+1) o ' o (—1), where ±1 have the obvious meanings. 
Thus 

H'W = fjfe'\ L; = L f and = A( M «, s^) = • s) Lfc - fc -i 
where k = (J — 1)'. 

Theorem 4.3 (Cf. JE/) Let n > 1 and s G D n and assume the initial values of Definition 
4-2\ Define jjS^ recursively by 



[ ti (n - 1] */A n = 

A* W = 

( A' n ■ x max K-i-o} ^(n-i) _ An . a jnax{-e B _i,o} ^/(n-i) otherwise. 

If A n = 0, clearly /J- n > G MP(s), L n = L n _i and e n = e n _i + 1. // A n 7^ taen 
(%) deg(/i (n) ) = max{e„_i, 0} + L n _i = n' + 1 — 14 

(n; g mp(s) 

(m) e n = -|e n _i| + I. 
Corollary 4.4 (Cf. JE/) For any s G D n , L n = n' + 1 - 14. 

Corollary 4.5 (Iterative Form of yu) Let n > 1 ; s G -D" and e £ D. Assume the 
initial values of Definition \4-S\ For 1 < j < n, let 



pV-V «/A i = 

A'j ■ _ ^. . x m^{- ej -xfi}^{3-i) otherwise. 



Then G MP(s (i) ). Further, if Aj = 0, then = /i'V-V, A' j+1 = A' 5 and e 3 - 
Cj-i + 1. If Aj 7^ inen 



(aj i/e^i < then = /Z^ 1 ) and A£ +1 = A$ 
f&) but ifej-i > f/ien = and A£ +1 = Aj 

(c) e,- = -|%_i| + 1. 



7 



In other words, when Aj ^ 0, 

f A' • fi^-V - Aj ■ x- e fi f( i-V if e = e,-_i < 

( Aj ■ x +e ^~^ - Aj ■ otherwise. 

See [10] for the derivation of the following algorithm from Corollary 14.51 

Algorithm 4.6 (Iterative minimal polynomial) 

Input: n > 1, e G D and s = (s±, . . . , s n ) G D n . 
Output: (i G MP(s). 

{e := 1; // := e; A' := 1; /i := 1; 
FOR j = 1 TD n 

IF A ^ THEN {IF e < THEN fi := A' • /i - A • x" e /i'/ 

ELSE {£ := //; 

/i := A' • x e // — A • 
^ ;= t ; A' := A; 
e := -e}} 

e := e + 1} 
RETURN (/i)} 

5 Minimal Realisations for D n 

We recall the definition of a minimal realisation from [12]; when /z G MP(s) is monic, this 
amounts to rational approximation of s by a rational function /^//-t for a certain /i 2 G 
with deg(/i 2 ) < deg(/x) and deg(/x) minimal. 

Recall that v : -D[x _1 ,x] — )■ Z is v(g) = m&x{i : ^ 0}. Let / G -D[x] x and 
d = deg(/). Then for s G D n , v(s) < —1 and 

v(s)+rf / d-n-1 -1 v(s)+d\ 

/•a= E(/-fiW= E + E + E 

j=-n yi=-n j=d-n j=0 J 



Definition 5.1 For fixed s G D n and f G D[x], / 2 = /2(s) G Z)[:r] is i/ie -Djxj-sttmmand 
off-s: 

v(s)+deg(/) 

/ 2 (x)= e (/-aW 

i=0 

and7=(/,/ 2 )G J D[x] 2 . 
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Thus / • s = (/ • s- f 2 ) + h, and deg(/ 2 ) = v(s) + deg(/) < deg(/). If / is monic, 
v(/)v(/ a //) = v(/ 2 ) and v(/ 2 //) = v(/ 2 ) - v(/) < i.e. f G af^aT 1 ]]. The next 
result is immediate, by construction. 

Proposition 5.2 7/s G -D n ; / G Ann(s) if and only if f = orv(f-s — / 2 ) < deg(/) — n. 

Thus if/G Ann(s) is monic, v(s-f 2 /f) = v(/-s-/ 2 )-v(/) < v(/-s-/ 2 ) <deg(/)-n-l 
and 

s-^ = 0mod (x^^-™" 1 ). 

Definition 5.3 (Minimal Realisation) f|7]| Sections 2,3]) We will say that f G D[x] 2 
realises s if either (i) f = or (mJ / 7^ 0, deg(/ 2 ) < deg(/) andv(f ■ s — / 2 ) < deg(/)— n. 
Further, f is a minimal realisation o/s i// realises s and f G MP(s) ; written f G MR(s). 

Then if / ^ 0, / G MR(a) if and only if / G MP(s). 

For //W as in Theorem 14.31 we write /4 f° r the polynomial part of //^ • s 

v(s)+d 
fc=0 

By definition, if yujf 7^ 0, then deg(//2 ) — v( g^l ) + deg(//^) < deg(/i^) and it is easy to 
check that deg(/i 2 J ^) < deg(/4^), where 1 < j ' < n. If A n = 0, then we obviously have 

(n) (n-l) 

It is convenient to separate out the following special case before extending Theorem 
14.31 to minimal realisations. 

Proposition 5.4 Let n > 1, s G -D n , L n _i = and A n 7^ 0. // /4 = — 1 anc ^ /4 = 0? 

Proof. We have to show that the second components of both sides agree. It is clear 
that s has precisely n — 1 leading zeroes, which gives the following data: Aj = for 
1 < j < n — 1, yifj^ = 1 for < j < n — 1 and /4" ^ = 0, e n _i = n, (n — 1)' = —1 and 
A n = s„, A^ = 1. Substituting these gives /4 = s n = A n ■ x n jjL^^ — A n • /4 , as 
claimed. ■ 

An iterative version of the following result (without the value of deg(/4"^)) was proved in 
|12j . [T4] . For completeness, we prove the recursive form. 

Theorem 5.5 (Recursive MR) (Cf. JIM/ ) Let n > 1 and s G -D n . Assume the initial 
values of Definition \4-2\ and set /i 2 X ^ = —1 and /4 = 0. Define recursively by 

( £ (n-1) */A n = 

[ • x max K-i.°} _ A n • xu**{-°»r-u0} p/(n-i) otherwise. 

If A n = 0, clearly jft 1 ' G MR(s), L n = L„_i and e n = e n _i + 1. If A n 7^ then 
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(i) deg(//W) = max{e n _i, 0} + L n _i = n' + 1 - L' n 

and z/ // 2 7^ i/ien deg(/4™^) = max{e n _i, 0} + deg^™^) 

(ii) /Z (n) e MR(s) 

e„ = -|e n _i| + 1. 



Proof. Put u = v(s), n = L = L„_ a , e = e n _ 1; // = /i'^ 1 ), L' = L^, 

A' = A' and A = A„. From Theorem 14.31 we have 



^ {n) = A' ■ x max{e ' 0} fi - A ■ x max {- e -°} y/. 
First we show that /i 2 (n) = A' • x max ^ /j, 2 - A ■ x max {- e >o} // 2 . For e < 

= ^ (^(») • s) 3 ;X j = A' ■ fl 2 - A • V • s)j- X J ' 

j=0 j=0 

since L n = L. The second summation is Y^k^e^ip' ' xk ~ e i- e - 



-1 v+L' v+h+e s 



v+L'+l/ 



say. Now for n > 1, (n-1)' < n-2 < n-2 + L'. From Theorem i ), L = (n— 1)'+1— L' 
which gives L' — (n — 1)' < e since L < n — 1 and so A = 0. Likewise, L + e = n — L = 
n - 1 - (n - 1)' + V > L' + 1 and so 5 = 0. Thus /4" } = A' • /i 2 - A • x~ e /x' 2 if e < 0. 
For e < 0, we have 



.(») 



D+n— L 



/i 2 ' = 2J ' a;e / i " £ — A • // • s) 3 - x- 7 

D+n— L n+n— L 

a 7 ■ ^ (x e /x • s)j- x 3 - a ■ 53 y • s 

i=o i=o 

(-1 u+L\ u+L' 
fe=-e fc=0/ j=0 
-1 

A' • x e E (A* ■ s)* a; fc + A' • x e fi 2 - A • // 2 

k=—e 



3 xJ 



since t> + n — L — e = v + ~L and u + n — L>u + L>u + L/. From Proposition 15.41 we 
can assume that L > 1. This implies that L — (n — 1) < — e and so the first summand 
vanishes, giving /i 2 n ^ = A' • x e ^ 2 — A ■ /i 2 if e > 0. 
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Next we prove that if /4 7^ then deg(/4 ) = max{e, 0} + deg^)- Since ^ 0, 
s (n-i) ^ (o, . . . , 0) and so v' = v(s n _ix 1 ~ n + • • • + siar 1 ) = v. Thus deg(/i 2 n) ) = n + L n = 
t>' + max{L, n — L} = t>' + L = deg(/Z2) if e < 0. If e > 0, then L n = n — L = L + e and 
deg(/4 n) ) = v + (L + e) = (7/ + L) + e = deg(/x 2 ) + e. ■ 



Corollary 5.6 (Iterative Form of //) Lei n > 1, s G -D n and e G D. Assume the 
initial values of Definition \4 . 2\ For 1 < j < n, let 

■pV- 1 '* if A d = 



A'j ■ x max { e .-i.o} ^j(i-i) _ . zMBxi-ej-ifl} j['U~i) otherwise. 



Then /Z^O G MR(s (j) ). Further, if Aj = 0, then jfi) = Jjfi-V ', A' j+1 = A'j and = 
e 3 -_i + 1. IfAj^O then 

(a) i/e^i < t/ien /I'Cj) = and A' j+1 = A'j 

(b) but if e^x > then Ji'^ = /Z^ 1 ) and A' j+1 = Aj 

(c) ej = -| e i-il + 1. 

For future use, we recall: 

Lemma 5.7 Lei s G D n . 

(i) ( ITfy, Proposition 2.3], \TJ\ Lemma 6.1]) If g G /i G Ann(s) and deg(g) + 
deg(/i) < n then {gh) 2 = gh2 . 

(ii) J/deg(/) < ra — L and deg(f') < L — L' then 

L'-l 

(/■/i + / , -/i / ) 2 = /^ 2 + / / - E (/*'■*)* x J - 

3=0 

Proof, (i) Let H = h-s — h 2 , d = deg(g) and e = deg(h). Since h G Ann(s), hs = H + h 2 
and (gh) ■ s = g(h ■ s) = gH + gh 2 . As v(gH) = deg(g) + v(#) <d-n + e- l<-l, 

Ma = ^2- 

(ii) Let d = deg(/). Then from Part (i), 

d+L d+L 

(f-v+f- = (/ • ^2 + 5^(/y • x j = / • /i 2 + j^(/v • s)j ^ 

j=0 i=o 
and so it suffices to evaluate the second summand. We have 

t-n+L'-l -1 L'-l\ 

E + E + E o'-sW 
j=—n j=— n+L' j=0 / 

Now deg(/') — n + L' — 1 < L — n — 1 < —1 since L < n and the middle summand is zero. 
Finally, if L' < j < d + L then (// ■ s)j = since deg(/x') = L' and s -_ L , =0. ■ 
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Proposition 5.8 ( [TR Proposition 4-13]) Let n > 1, s G D n , f G D[x] and deg(f') < 
—e n . Then jtf- n ' + /' ■ Jl'^ is a minimal realisation of s. 

Proof. For a proof that fj,^ + f ■ fx'^ G MP(s) using the index function, see [HI 
Proposition 4.9]. The proof that we have a minimal realisation now follows from Lemma 

o ■ 

Theorem 5.9 (Jig Theorem 4. 16]) Let n>l, D bea field and s G D n . IfJ G MR(s), 
then f = Jl^ + f ■ ~pf( n > for some f G D[x], where f = or deg(f') < —e n . In particular, 
if e n > 0, then /Z^ is unique. 

5.1 The Iterative Algorithm 

As for the minimal polynomial algorithm, we write 7Z'^ for jft' and then suppress j. 
Algorithm 5.10 

Input: n > 1 and s = (si, . . . , s n ) G -D n . 
Output: ~fi G MR(s) and rfs Mi? identity. 

{A':=1/7Z':= (e, -1); e := 1/ /Z := (1, 0); 
for j = 1 TO n 

IF A ^ THEN {IF e < THEN IZ := A' • 7Z - A • x" e 7Z' 

ELSE {t := IZ; 

/Z := A' • x e /Z - A • /Z'; 
7/ : = t; A' := A; 
e := -e}} 

e := e + 1;} 
return (/Z)} 

Example 5.11 (Cf. /Iffi Table I]) Table \5.1\ gives the values ofjl^ andJjfV'. The reader 
may check that p/V' ■ jfo) = 1 for 1 < j < 8 as per Theorem \6.8l 

Over a field, we can produce a monic \i by dividing by A'. We remark that \x and ^2 can 
even be computed in parallel. 
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Table 1: Algorithm EM with input (0, 1, 1, 0, 0, 1, 0, 1) G F§. 



3 


A? 


e i-l 








Mil 





i 




1 








1 


1 





1 


1 








1 


2 


i 


2 


x 2 


1 


1 





3 


i 


-1 


X 2 + X 


1 


1 





4 


i 





x 2 + X + 1 


1 


1 





5 


i 


1 


X 3 + x 2 + X + 1 


x 2 + X + 1 


X 


1 


6 








X 3 + x 2 + X + 1 


x 2 + X + 1 


X 


1 


7 


1 


1 


x 4 + X 3 + 1 


X 3 + x 2 + X + 1 


x 2 + 1 


X 


8 


1 





4 i 2 i 
X + X + X 


X 3 + x 2 + X + 1 


x 2 + X + 1 


X 



6 Minimal Polynomial Identities 

6.1 An Identity for /i^ and 

Our first set of identities was suggested by the following result. 

Proposition 6.1 Let D be a factorial domain and s G D n . For 1 < j < n, let e = Qj-i. 
Then 

(i) g<A<jiV\n'®) = \ 

(ii) if Aj 7^ and pj > is the highest power of x dividing fi^\ then 



1 i/e>0. 



fzzzj gcd^W*,^'-^*) = 1. 

Proof. We prove this by induction on j, the case j = being trivial. Suppose inductively 
that gcdOxG- 1 ),^'- 1 )) = 1. Then 



/' 

and 



(i) _ J A i ' - A i • x~ e ^' {j ~ 1] if e < 

if e < 



A' ■ x +e ^ j -V - Aj ■ ^'C*'- 1 ) otherwise 



// 



'(j) - J V 



if e > 0. 

Thus if e < 0, then gcd(/i^,//^) = gcd(>^, /Z^" 1 )) = gcd^' -1 ), p,'^) = 1 by the 
inductive hypothesis. If e > 0, then 

gcd(/z w) ,// a) ) = gcd{fj^\ = gcd^'G- 15 ,^- 1 )) = 1 

by the inductive hypothesis. 
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(ii) Suppose that e < 0. Then fi® = A' f fj,^-A r x- e /j,'^ and so gcd^^ 1 ), /ft)) = 
gcd(/i^- 1 ), x~ e n'^) = xriato-u-e} by Part (i). For e > 0, we have (i® = A'j ■ x e /x^ - 
Aj ■ n'v~ l > and so 

gcdO '- 1 ),//^) = gcd(^ j - l) , = 1. 

(iii) Easy consequence of (ii). ■ 



Corollary 6.2 Let D be factorial and s G D n . The successive 'connection' polynomials of 
s (obtained by the Berlekamp-Massey algorithm for D of FfU, Algorithm 5.4]) are relatively 
prime. 

Proof. For 1 < j < n, the successive connection polynomials of s are [$) * and jtfk'- 1 ) * 
by [TOi Corollary 4.10]. The result now follows from Proposition 16. l( iii) . ■ 

We will see that the following inductively defined scalar plays a key role: it turns out to 
be the right-hand side of all three identities. 

Definition 6.3 Let s G D n and for 1 < j < n, let Aj, and e = e^-i be as in Theorem 
15.51 Put Vo = 1 and for 1 < j < n, define V j = Vj(s) by 
(i) V^Vi-i ifA s = 

(ii) 

A', ■ Vj_i ife<0, A^O 



■ Vj-i i/e > 0, Aj ^ 0. 



We have Vi = Ai since eo = 1. It is easy to see that the Vj are always non-zero. We 
now construct an identity for /r n ) and (j,'( n \ 

Theorem 6.4 Let D be a domain, s G D n and for 1 < j < n, let e = ej_i be as 

(0) (fc) 

in Theorem Put f = (1,0). Assume inductively that we have defined f for 
< k < j - 1. IfAj = A.-^O-D) = o, put f j) = 7 0_1) . // Aj + 0, pit* 



7 0) 



(+/(^- 1 ), A;. ■ /i^ + Aj ■ x-^/C?- 1 )) i/e < 

(-/ 2 °'- 1) ) a;. • ,••'./;? 1 + Aj • i/ e > o. 

JTien /or < j < n ; we /iawe /^ ■ (y/^,//^) = Vj. 

Proof. The case n = is a trivial verification. Suppose that /'' • (p^\ /j,'^) = 1 for 

< j < n - 1. Put /i = ^ n ~ l \ e = e n _ 1; A = A„, /V = ^ n ~ l \ A' = A' n and / = J (n_1) . 
We can assume that A„ 7^ 0. If e < 0, we have = fj/ and 

/ (n V n) = /(A' • /x - A • x- V) = A' • f/x - A ■ x- e fix' 
= A' • (Vn-i - /2/i') - A • x-'f/j,' 
= V n -(A'-f 2 + A- x- e f)fi' = V n - / 2 (n V (n) - 
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If e > 0, we have /i /(n ) = [L. Thus 

/ (n V n) = -/ 2 (A'-x e /i-A- /U ') = -A / ^72/i + A-/ 2/ i' 

= -A' • x e f 2t i + A ■ (V„_! - ffi) = -(A' ■ x e f 2 + A ■ f)n + V 

= -/jV B) + V„. 



When £) is a field, these amount to Bezout Identities for the fi^ and //( n ). We see again 
that gcd = 1 if D is factorial. The general form of /g shows that we can 
compute it iteratively alongside in Algorithm 15. 101 Our first application follows from 
Proposition 16.11 when D is factorial. 

Corollary 6.5 Let D be any domain and s G D n . If //j, = then //q ^ and similarly 
for fio 

If needed, we can obtain deg(/ 2 (n) ) from deg(/( n )) + L n = deg(/ 2 (n) ) + L' n . 
6.2 An Identity for ^ n_1 ) and ^ 

We give an identity for // n_1 ) and fj,^ when A n ^ and e n _i > 0. When A n 7^ 0, 
e n-i < and D is factorial, there is a similar non-constructive identity (which does not 
explicitly involve V„). 

Proposition 6.6 Let s G D n , A n 7^ and e = e n _i be as in Theorem \5.5\ If e > and 

J (n ' 1] ■ (/i^- 1 ),//™- 1 )) = V n _i ^en 

(A ■ /C"- 1 ) + A' • x°ft l \ ~ft l) ) ■ ^ n - 1] ^ {n) ) = V„. 

Proof. Put fi = n^-V, e = e n _i, A = A n , // = ^ n - l \ A' = A' n and J = J M '. The 
identity is trivial if n = 1. Suppose that f^ 1 • (^ n ~ 1 \ fi r ^ n ~ 1 ^) = 1. Then 

(A./ + A'-x e / 2 ,-/ 2 )-(^^ (n) ) = A-//i + A'-a;72/i-/ 2 /i (n) 

= A • (V n _! - / 2 //) + A' • x e f 2l x - f 2 »W 

= V n — A • / 2/ / + A' • / 2 x e /i - / 2 //W 

= V n + f 2 (A'-x e fi-A- l 2 , )-f 2 ^ = V n . 



The case e„_i < is complicated by gcd(/i^ n e ) = x m say. If m = — e, one may 
check that 

In general, we can apply the following result when D is factorial. 
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Theorem 6.7 (JEi Corollary 2.7]) Let D be factorial and a,b G D[x], not both zero. 
There are u G D and f,gE D[x] such that fa + gb = ngcd(a, b). 

For a specific s G D n , we may even obtain u, f and g by applying algorithm XPRS of 
[3] to /Z"' -1 ) and See also Theorem 18.11 It would be useful to have an identity for 
general D which involves V n when A n ^ and e n _i < 0. 

6.3 An Identity for /x( n ) and fi^ 

Recall that Jl^ = (e, -1), Jl^ = (1, 0) and that if / = (/, / 2 ), then / = (-/ 2 , /). 

Theorem 6.8 For < j < n, we have ft'^ - ~p^ = Vj. Hence if ^ = then (/4 )o 7^ 
and similarly for (/4^)o- 

Proof. For j = 0, we have ■ (1,0) = 1 = Vo- Suppose 1 < j < n and the result is 
true for k < j — 1. If Aj = 0, there is nothing to prove. Otherwise, let e = Qj-i- If e < 
then 

= -/i' 2 (A;. ■ [i — Aj ■ + //(a; . M2 - a, • x-y 2 ) 
= A' r pf--p = a; ■ v,-_! = v, 

and if e > then 

= -/i 2 (A^ • x e [i - Aj ■ (j!) + n(A'j ■ x e ^ 2 - Aj ■ fj! 2 ) 
= Aj ■ x e {-\i 2 [i + /i/i 2 ) + Aj ■ (-fjhfJL + ///i 2 ) 
= Aj ■ ft -~p= Aj ■ Vj-i = Vj 

which completes the induction. ■ 

When D is a field, the identity for /i 1 -™ -1 ) and n^ n ' of Theorem 16.81 amounts to a Bezout 
identity. The following result was proved differently in jX2|, Corollary 3.24]. 

Corollary 6.9 Let D be a factorial domain and s G D n . Then gcd(/7^') G D for 1 < 
j < n. 

Corollary 6.10 Let D be a field and s G D n . Then for < j < n, p,'^> are Bezout 
coefficients of any minimal realization of s. 

Proof. Let m be a minimal realization of s. By Theorem I5.9[ m = ~p,( n > + m'jl'^ = 
Ji + m'~p' say, where deg(m') < — e n . Now fi! ■ ~p' — gives 

p! ■ m — fx' • ~p, + m'fjf ■~fi' = fi'-JI = V„. 
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Proposition 6.11 Let D be a field, s G D n and Bezout coefficients a be computed for 
/j( n ) by the extended Euclidean algorithm. Then V„ • a = Jl'^ 1 ' . 

Proof. Put a = (a, a 2 ), Jl = J^ n \ j! = ~p'^ and V = V n . We have a ■ Jl = gcd(TZ) = 1 = 
V -1 • ft! ■ Jl. Hence (a + V _1 /4)yU = ~( a 2 — V -1 //)/-^- If — V -1 /^ = 0, we are done. 
Otherwise /x|(a 2 — V -1 //) since gcd(/l) = 1. This is impossible since deg(a 2 — V~V) < 
deg(». Thus a = V"V. ■ 

We can also apply Theorem 16.81 to linear recurring sequences. 

Proposition 6.12 Let D be a principal ideal domain and s G D n be a linear recurring 
sequence with minimal polynomial \x. There are explicit f G -D[x] 2 and V G D\ {0} such 
that f ■ Ji — V. 

Proof. Apply Theorem EH to s (2deg(/i)) . ■ 
6.4 The Iterative Algorithm 

The following extension of Algorithm 15.101 is immediately implied by Theorem 16.81 
Algorithm 6.13 

Input: n > 1 and s = (s%, . . . , s n ) G D n . 
Output: Ji G MR(s) and the MR identities of s. 

{pf:= 0,-1); A' := 1; e := 1; Jl := (1,0); 

V := 1; 

FOR j = 1 TO n 

{A : = E fe =o^ s fc+^; 

IF A^O THEN {IF e < THEN {ji := A' • /Z - A • x" e /Z'; 

V := A- V} 
ELSE {£ := 77; 

71 := A' • x e Jl — A • 77'; 
77' : = t; A' : = A; 
V:= A-V; 
e := -e}} 

e := e + 1;} 
RETURN (77,/i', V);} 

We could also compute Vj at the end of iteration as / • ~pP\ but this would increase the 
complexity of the algorithm. 
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Table 2: Theorem El with s = (1,0,1,1,0,1) G GF(2) 6 . 



3 


e 


A 




71' 




1 


1 


1 


(x,l) 


(1,0) 


(0,1) 


2 








(x,l) 


(1,0) 


(0,1) 


3 


1 


1 


(x 2 + 1, x) 


(x,l) 


(l,x) 


4 





1 


(x 2 + X + 1, X + 1) 


(x,l) 


(l,x) 


5 


1 





(x 2 + X + l,x + 1) 


(x,l) 


(l,x) 


6 


2 





(x 2 + X + 1, X + 1) 


(x,l) 


(l,x) 



Example 6.14 Let s = (1, 0, 1, 1, 0, 1) G GF(2) 6 . Then 

7 (3) • 7I (3) = 1 • (x 2 + 1) + x • x = 1 
7 0) . jft) = 1 • (x 2 + x + 1) + x • (x + 1) = 1 for 4 < j < 6. 

7 Perfect Linear- Complexity Profile 

The following generalisation of a definition from [16J to any domain D makes sense by 
Theorem 14.31 



7.1 Basics 

Definition 7.1 A sequence s G D n has a perfect linear-complexity profile (PLCP) if 
Li =[^r\ forl<]<n. 

It is easy to see that for D = ¥2, the sequences of length 1 to 4 with a PLCP are (1), 
(1, S2), (1, 1, 0), (1, 0, 1) and (1, 1, 0, S4), (1, 0, 1, S4). See Table 17.11 2 for the corresponding 
~pV\ where Aj = A( / u^ _1 ^). Recall that for any sequence, p ' = 1 and 6q = 1. 

Proposition 7.2 The following are equivalent: 
(1) s has a PLCP 
(ii) for 1 < j < n 

1 if j is even 



otherwise, 

(Hi) Aj 7^ for all odd j, 1 < j < n 
(iv) = (x + e, 1) and for 2 < j < n 



Aj_i • fj,V x ) — Aj ■ 1$ 2 ^ if j is even 



Aj_ 2 ■ x/i^ ^ — Aj ■ [ft 3 ^ otherwise. 
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Proof, (i) (ii): Easy consequence of the definitions. 

(i) =>- (iii): If j < n + 1 is odd then Aj ^ 0, for otherwise ^ + 1 = ^ = Lj = Lj_! = 

3-1 
2 - 

(iii) (i): Let Aj 7^ for all odd j, 1 < j < n + 1. Then Sl 7^ 0, L x = 1 and 
ei = 0. If A2 = 0, then L2 = Li = 1, otherwise L2 = max{ei, 0} + 1 = 1, so that L2 is as 
required. Suppose that j < n is odd and L& = [^"J f° r all /c, 1 < /c < j — 1. We have 
Lj = j— Lj_i = j — ^- = L^irJ • If J = ""-+1, we are done. Otherwise, if Aj + i = 0, we have 
L j+1 = Lj = L^J = L^J, whereas if A j+1 ^ 0, L i+1 = j + l-L, = j + L^J = L^J • 
(ii) =>■ (iv): We have 

= A' • x max { c J-i'°} — Aj ■ x max {- c 3-i. } -p'(i~ 1 ') 

where Aj may be zero and 6j_i is as usual. If j is even then Aj_i 7^ and ej_2 = 1, so 
i = j — 2; if j is odd then ej_ 2 = and ej_ 3 = 1 so (j — 2)' = j — 3. Inserting e 3 -_i now 
gives the formulae for ~p. (Note that (j — 1)' + 1 is always odd, so that A^- 7^ 0.) 

(iv) =>• (ii): We have L x = 1, Lj = Lj_i if j is even and Lj = Lj_i + 1 if i is odd. Thus 
if j is odd, ej = j + 1 — 2Lj = j + 1 — 2(Lj_i + 1) = j + 1 — 2Lj_i = e 3 -_i + 1. Applying 
this inductively gives (ii). ■ 



Remarks 7.3 (%) We applied Theorem \5.5\ to s G -D n , no£ to (Ai, . . . , A n ) G -D n as stated 
in JEQ. 

(ii) If s is an infinite linear recurring sequence [i.e. s G is a rational function 
gj f with deg(g) < deg(/)/ then lim^ooen = 00 whereas if s has a PLCP, trivially 
lim n ^ c ^e„ = ; lim n _>. 00 e n = 1 but limn^ooen does not exist. So linear recurring sequences 
never have PLCP. 

(iii) For an infinite sequence s, Proposition \ 7. 2\ is consistent with ' t 9l. Let 

K(s) = sup{deg(A n )} 

n>l 

as in [9]. Then if s is not a rational function and has partial quotients A i; Theorem 1 of 
/2J/ implies that for all n, 1 — K(s) < e n < K(s) and s has PLCP if and only if s is not a 
rational function and K(s) = 1/3 Theorem 3]. 

Proposition 7.4 If D — ¥ q , then \{s G D n : s has PLCP}\ = (q - l^flgLfJ, and if 
\D\ = 00 then {s G D" : s has PLCP} is infinite. 

Proof. A sequence s determines (Ai, . . . ,A n ) G D n uniquely, and conversely. Thus the 
result follows from Proposition I7.2( iii). ■ 

It follows that without loss of generality, we may assume that n is odd. It is easy to see 
that for odd n, /I (n) + c • /Z /(n) = /7 (n) + c • /Z^ 1 ) G MR(s) for any c G D. 
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Table 3: a® G ¥ 2 [x] for < i < 4. 



J 









(1,0) 


1 


1 


(x,l) 


X + 1 


2 


(x + A l5 l) 


(A x + l)x + l 


3 


(x 2 + A]X + 1, x) 


(Ai + l)x 3 + x + l 


4 


(x 2 + (Ax + A 3 )x + l,x + A3) 


(A 1 + l)x 3 + (A 1 A 3 + l)x + l 



Corollary 7.5 Let D be a field and fj,W and Aj = A^/i^ ^) be as in Theorem 15. .91 for 
l<j<n. If s has PLCP, then 



1 {fj,ti ^ — ■ fj,ti 2 **} if j is even 

MR(s (j) ) = I 

I {(x + c) ■ /7^ _1 ) — -^-^ ■ /l( J " 3 ) : c G -D} otherwise. 



Proof. If j is even, then s^ has a unique monic minimal realisation by Theorem 15.91 
Similarly if j is odd, any monic minimal realisation of s^' is 



for some c G -D. 



Corollary 7.6 7/s /jas PLCP then for 2 < j < n, then gv> ■ ^C?) = y . where 

{—^2 2 \ t 1 ^^) if 3 i s even 
(— , jU^' -1 )) otherwise 



and 

(Aj_i ■ Aj_ 3 • • ■ A1) 2 z/ j zs even 



A,- • (Aj_ 2 • Aj_ 4 • ■ ■ Ax) 2 otherwise. 



Proof. We have gU> = jx'^' and j' = j — 2 if j is even and j' = j — 1 if j is odd. To 
verify the values for Vj, we have Cj = 1 if j is even and Cj = if j is odd. Thus if j > 2 
is even, Vj = Aj_i • Vj_i and if j is odd, Vj = Aj • V 3 -_i, which easily yields the stated 
formulae. ■ 

With D = ¥ 2 for example, (0, 1) ■ (x + A 1 , 1) = 1 and (1, x + A x ) ■ [x 1 + A x x + l,x) = 
(x 2 + Aix + 1) + x(x + Ai) = 1. 
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7.2 Binary Sequences 

Here we give a simpler proof of a theorem of Wang and Massey (20] on binary sequences 
via Bezout coefficients. Thus D = F 2 throughout this subsection. Let us call s G D n 
stable if Si = 1 and for even j, 2 < j < n, Sj + \ = Sj + S£. It is convenient to use two 
auxiliary functions 

t = t(s) = s 2 + (x + l)s + 1 

and 

a (j) = u m + ( x + + 

where /7^) = (//to), i/Wj j s as j n Proposition I7.2( iv) ; t is also used in the continued fraction 
treatment of this theorem in [H] (for infinite sequences) . First we have the following simple 
consequence of the definitions. 

Proposition 7.7 fj,( n > 2 t = + F where v(F) < if n is even and v(F) < 1 if n is 
odd. 

Lemma 7.8 Let n be odd. The following are equivalent: 

(i) s is stable 

(ii) tj = for j even, < j < n. 

Proof. We have t = s\ + l and for all even j, 2 < j < n — 1, we have tj = Sj + Sj + \ + s i . 

■ 

The innocuous-looking part (ii) of the following lemma is essential to our proof. 
Lemma 7.9 (Cf. J2(% Lemma lj) If s has PLCP then (i) for 2 < j < n, 

{ a (j-i) _|_ A; • a( j ~V + A,- • (x + 1) if j is even 
x 2 a^~^ + cr^ -3 ) + x(x + 1) otherwise 

(ii) for < j < n, = 1 and = 0. 

Proof. Put v = /i 2 . If j is odd, ~pW = xjl^' 1 ^ + Jt^ j ' 3 \ so 

a (i) = [ x 2 zy ( J -l)2 + ^-3)2] + ( z + !) + M 0-3)] + 

+ [ x 2 /i (i-l)2 +/i ( J -3)2 ] 

We have /Z^ 1 ) = ZZ (j - 2) + Aj_i • 7Z (j '~ 3) , so 

= V j_3) + ^-V^'- 3 ) + 2A i _ 1 • vV-QpV-Q 

= ^-2) M 0-3) + ^0-2)^-3) = x 
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by Corollary 17.61 and so for i odd, = x 2 a^ ^ + 3 ' + x(x + 1). If i is even, 
■pfi) = jfi-i) + A . . jfi-q and Corollary E6] yields 

' +[^ 2 + A l -^ 2 ] 

= + Ai • a^-V + Ai ■ {x + 1) 

by Corollary 17.61 We omit the simple inductive proof of Part (ii). ■ 

The reader may also check that deg(c^) = j — 1 if j is even and j if deg(cr^) = j is odd, 
but we will not need this fact. 

Theorem 7.10 Let s G D n . Then s has PLCP if and only if s is stable. 

Proof. It suffices to show that if n is odd and s has PLCP then tj = for even j, 
< j ' < n — 1. Put // = //W). Inductively, we have 

/i 2 • (tix" 1 + • • • + t n - 2 x 2 ~ n + t n _ x x x - n + t n x~ n ) = /j 2 -t = a + F 

where v(F) < 1. We want to show that t n -\ = 0. The coefficient of x 2 in the left-hand 
side is /in+i ■ t n _i and a 2 = by Lemma ITTW ii). Since is the leading coefficient of 
/i, we must have t n -\ = . Conversely, as in noted in [20], there are clearly |~~] stable 
sequences in Z) n , so the result follows from Proposition 17.41 ■ 

8 An Algorithm for the Polynomial Bezout Identity. 

Theorem 8.1 (Cf. /3J Corollary 2.7]) Let D be a principal ideal domain, u G -Dja;] 2 , 
with u monic and < deg(«2) < deg(w). There are explicit f G D[x] 2 and V G D \ {0} 
such that f -u = V ■ gcd(w) . 

Proof. If deg(w 2 ) = deg(w) = d, let £, £ 2 be the leading coefficients of u and u 2 respec- 
tively. We replace u 2 by u' 2 = 1% • u — I • u 2 . If / • (u, u' 2 ) = V • gcd(w), one checks that 
(/ + ^2/2, —£f2) • u — V • gcd(w). So we can assume that deg(u 2 ) < deg(w). Let t be the 
linear recurring sequence defined by t = I G D[[ X -% s = t^ and /Z be the unique MR 
of s, with /, V as in Proposition 16.121 so that //Z = V. 

We know that u G Ann(s) = fiD[x] and soti = wfi for some w G D[:r]. We show that 
w = gcd(w). For u 2 = ut = {wjj)t = w(fit) = wfi 2 i.e. u — w ■ ~p, so gcd(U)|u\ We also 
know that there exists a G Dfx] such that a • u — gcd(w). Then gcd(w) = a ■ w~p, = wa ■ Ji, 
so w \ gcd(w). Finally, V • gcd(w) = V • w = w ■ (/ • /I) = / • (w~p) = fu. ■ 

As in Proposition 16. 11[ a degree argument shows that V -1 / agrees with the coefficients 
found by the extended Euclidean algorithm. Thus Algorithm 18. 21 will be widely applicable, 
not just to ¥[x], ¥ a field, but also to ¥[x, y] for example, as ¥[x] is a principal ideal domain. 
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Algorithm 8.2 (Cf. [3, Algorithm 4.2]) 

Input: P.I.D. D, u G -D[x] 2 , with u monic and < deg(w 2 ) < deg(tt) = d. 
Output: f G -Dfx] 2 such that f ■ u = V • gcd(u). 

1. Compute s = f( 2d > by 2d subtractions of u from U2 in D[x~ l ,x\, where t=—. 

2. Apply Algorithm \6.13\ to s, giving f G D[x] 2 and V G D. 

3. Compute f ■ u, which is V • gcd(w). 

Note that step 2 requires at most 5d 2 multiplications (and at most 3d 2 multiplications if 
D is a field). 

Example 8.3 Let D = GF(2), u = (x 3 + l,x 2 + 1) with gcd(w) = x + 1. Then 2d 
subtractions ofufromu2 in D[x _1 ,x] gives the generating function x^ 1 + x~ 3 + x~ 4 + x~ 6 
of s i.e. s = (1,0,1,1,0,1) as in Example \6.14\ We get 1 • (x 3 + 1) + x ■ (x 2 + 1) = x + 1. 



9 Annihilators Which Do Not Vanish At Zero 

First we establish a lower bound lemma which was stated without proof in |18J . The 
resulting easy corollary and algorithm do not require any characterisations of minimal 
polynomials. Theorem 16.51 (or Proposition 16.11 if D is factorial) allows us to remove a test 
on \j! . We also generalise Proposition 15.81 and Theorem 16. 8l to give another construction of 
these annihilators (by extending the original sequence by one term) and include minimal 
polynomial identities. In the final subsection, we prove a characterisation which was also 
stated without proof in [18] . 



9.1 The Lower Bound Lemma 

First of all, s G D n and we can assume that s is not the all-zero sequence, so that L(s) > 1. 
It is convenient to write Ann*(s) = {/' G Ann(s) : fo ^ 0}. We seek an / G Ann*(s) such 
that deg(/) = min{deg(g) : g G Ann*(s)}. To simplify the notation, we write JI = ~pS n \ 
e = e„ > and L = L n , each obtained as in Theorem 15.51 We also write L' = L n /, e' = e n >, 
Jl' = Ji {n '\ 

MR'(s) = {(/, f 2 ) : / G Ann'(s), deg(/) is minimal}, 

L* = L*(s) = min{deg(/) : / G Ann*(s)} and for n' > 1, s' = (s 1; . . . , s n /). 

If fi 0, clearly ~p G MR'(s) and evidently MR'(si, . . . , s n+1 ) C MR*(si, . . . , s n ). 
The case e < is an easy consequence of Proposition 15.81 but the case e > requires 
some preparation. We begin with the following easy result. 

Proposition 9.1 Let f G Ann(s) and d = deg(f). 

(i) Ifl<k<d and x k \f then f /x k G Ann(sfc + i, ...,s n ). 

(ii) If 2 < k < n and g G Ann(sfc, . . . , s n ), then for any t G D k ~ l , x k ~ 1 g G 
Ann(ti, . . s*., . . . , s n ). 
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Proof, (i) Let /' = f fx and s' = S2X~~ l + ■ ■ • + s„x 1_n . For d — n — (d — 1) — (n — 1) < 
j < — 1, (f's')j = (fs)j — S\fj + i = since / G Ann(s) and j + 1 < 0. Now induct on k. 

(ii) If g G Ann(s) then for any t G D , xg G Ann(t, si, . . . , s n ): for d + 1 — (n + 1) = 
d — n < j < —1, (xgitx -1 + Six~ 2 + ■ • • + Sn^ - ™ -1 )),,- = 9jt + (g§)j = since j < — 1 
implies that Qj = and g G Ann(s). The result now follows by induction on k. ■ 

Lemma 9.2 Let n > 2 and g G Ann*(s). 

(i) If h G Ann(s 2 , . . . , s n ) \ Ann(s), then deg(g) >n — deg(/i). 

^ // / G Ann(s), a;|/ and f/x£ Ann(s( n-1 }) ; i/ien deg(g) > n + 1 — deg(/). 

Proof, (i) Let d = deg(g). Then (gs)j = for d — n < j < — 1 and so we can write 
gs = G + ex" 1 -" 1 ' 1 + P where G G Ufar 1 ], Gj = for d - n - 1 < j < -1, c G £> 
and P G P[x]. Put s' = (s 2 , • • • , s n ) and e = deg(h). Then we have (fosQj = for 
e — n < j < —2 and a = (/is^)_i 7^ 0. As s = x~ V + six -1 , we get 

{0 if e - n - 1 < j < -3 

a if j = —2 

6 = (/wOo + si/ if J = -1 

and so hs = H + ax~ 2 + bx^ 1 + Q where H G Pfx -1 ], Hj = for e — n — 1 < j < — 1 and 
Q G -D[x]. This gives /iP — = gP — /iG + gax~ 2 + g&x -1 — hcx d ~ n ~ x G P[x] and 

= (#P)- 2 - (^G)_ 2 + #0^ + h n+ i^ d c. 

If d + e — n — 1 < —2, then n + 1 — d > e and so this reduces to = ago, for a contradiction. 
Hence d + e — n — 1 > — 1, as required. 

(ii) Proposition 19 . 1 If i) implies that h = f/x G Ann(s 2 , • • • , s n ) and as Ann(s) C 
Ann(s^ n_1 ^), h G" Ann(s). Part (i) now implies that deg(g) > n — deg(/i) = n + 1 — deg(/). 

■ 

The next corollary was stated without proof in jTH Proof of Theorem 3.7] when D is 
a field using different notation. 

Lemma 9.3 Let /x = and e > 0. If g G Ann'(s) i/ien deg(g) > n + 1 — L. 

Proof. We have n + 1 — 2L = e > 0. If s is the all-zero sequence, /x = 1 and //q 7^ 0, 
so s is not the all-zero sequence and 1 < L < |. Next we show that fi/x ^ Ann(s^ n ~ 1 ^). 
For suppose that fi/x G Ann(s( n_1 ^). Then L n _x < L — 1 = max{L n _ 1 ,n — L n _i} — 1, 
which implies that L n _i < L — 1 = n — L n _i — 1 i.e. e„_i > 0. Now Theorem 14.31 implies 
that e = — e n _i + 1 < 0, which is a contradiction. Thus fi/x G" Ann(s'" -1 ') and the result 
follows from Lemma l9.2( ii). ■ 
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9.2 Two Corollaries, an Algorithm and a Theorem 

Recall that e = e n = n + l — 2L. 

Corollary 9.4 (Cf. JM Theorem 3.8]) If fi = 0, let q G D[x], deg(g) = e, a e D x and 



Then (i) ft G MR*(s) and (ii) fx! ■ fx' — V n if e < 0; fx ■ fx* = —a ■ V n if e < 0. 

Proof. Corollary 16.51 implies that ft' ^ 0, so fx' ^ 0. If e < 0, Part (i) follows from 
Proposition 15.81 with f — a and Part (ii) follows from Theorem 16.81 Let e > 0. We check 
that fx' G Ann(s) directly from the definition. Let e + L — n — 1 — L < j < — 1. Then 

(A 4 * " s)j = (q ■ ^ ■ l)j + a- ■ <y • s) j = 

since (i) 1 — L < j implies that L — n < j — e and (q ■ p, ■ s)j = YH=o9k(p • §_)j-k, where 
j — e < j — k < j and (ii) n'—U < j. We have deg(/i*) = e+L, which is minimal by Lemma 
19.31 Finally, L > 1 implies that e < n — L and hence by Lemma [5.7( h) . = (?A*2 + a • ^2; 
so that ft* = q -~p + a -~p' G MR*(s). For Part (ii), we easily get jl-Jt* — —a ■ p! -~p which 
is — a ■ V n by Theorem 16.81 ■ 

A trivial example of the first case is n = 1 and S\ G D x . We have fx = (x, 1), fx' = (1, 0) 
and e = 0. According to the Corollary, we can take fx' — {x + a, 1) where a G D x . For 
the second case, let n — 4 and s = (0, 1,0,0). Then fx = (x 2 , 1), /I' = (1,0) and e = 1. 
According to the Corollary, /x 9 = x e fx + a(l, 0) = (a; 3 + a, x). 

Lemma [9.31 now yields the following value of LV Recall that if e < then fx G MP(s) 
is unique. 

Corollary 9.5 For s G D n , 



We now have all the ingredients for an iterative algorithm for sequences over D by taking 
q = x e and a = 1 in Corollary 19.41 

Algorithm 9.6 (Cf. JM Algorithm 3.2]) 

Input: integer n > 1 and s — (si, . . . , s n ) G D n . 
Output: fx' G MR'(s). 

{Algorithm \5.10\ (e = 0, n, s, fx, fx', e); 

IF /i /0 THEN /Z* := fx 

ELSE {IF e < THEN fx* :—fZ + fx' ELSE fx' := x e ■ fi + fx'}; 
RETURN (fx')} 




fx + a ■ fx' if e < 
q -Ji+ a - ft otherwise. 




n + 1 — L otherwise. 



if e < or (e > and /x 7^ 0) 
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Remarks 9.7 (Cf. fT8\ Section III]) (i) The correctness of Algorithm \9.6\ does not require 
us to characterise elements of Ann* (s) of minimal degree, (ii) Theorem \4-3\ uses any 
^(n-i). ^ n p ar ti cu l ar ^ if 1 < j < n, e 3 _x < and [Iq = ; we can replace ZZ^ -1 ) by 
jjU- 1 ) _|_ t/O'- 1 ) . (mj Algorithm \9.6\ does not include any tests on fi' . (iv) We can also 
use e = 1. 

Example 9.8 We apply Algorithm \9.b\ to Example \5.11\ ZZ® = (x 4 + x 2 + x, x 2 + x + I), 

-p/{8) _ ^3 _|_ x 2 _|_ x _|_ x ^ ana i 6g = SQ -p(8)m _ x jj(8) _|_ jj(S') = ^5 _|_ x _|_ ^ ^3 _|_ ^ _ 

The following theorem is our analogue of Proposition 15.81 and Theorem 16.81 
Theorem 9.9 (Cf. [18]) Let s G D" and f G D[x}. 

(i) If f = o or deg(f) < -e and /x + ft ■ /i' ^ then (a) JZ* = /Z + /' • // G MR'(s) 
and (b) fi' = V n . 

(ii) If e > and /i = ; /ei s n+ i 6e stzcn t/iat A n+1 = A„ +1 (/i,t) 7^ 0, where 
t = (si, . . . , s n , s n+ i) . If f = or deg(f') < e - 1 tnen 

(a) jr = /Z (n+1) + /' • /Z G MR*(i) C MR'(s) 

(b) ii-jr = V„+i(t). 

Proof. As in Proposition 19.41 we know that /Xq 7^ by Proposition 16.41 so the first 
part of (i) is clear from Proposition 15.81 and Theorem 16.81 and /i' • JZ* = fi' ■ ~fi- To 
prove Part (ii), first note that by Theorem |4J3l fi^ n+ ^ = ■ x e /z - A n+ i ■ (j? G MP(t) 
with deg( / u (rt+1) ) = n + 1 - L and A' n+1 = A n / +1 . Further, /i[, n+1) = A n+1 • fi' ^ 
and Ann(t) C Ann(s), so Lemma implies that ZZ (n+1) G MR(t) n MR'(s). We have 
e n+ i = -e + 1 < 0, so Zl( n+1 ) + /' • /Z' (n+1) g MR(t) by Proposition EHJ Also, (n + 1)' = n 
since e > 0. We conclude that = ^ 0, p* = ZZ (n+1) + /' • JZ G MR'(t) and 

/2 . /Z- = /x'^ 1 ) • ZZ* = • /Z (n+1) = V n+1 . 

■ 

In Example 19.81 S9 = forces A9 = 1. We obtain ZZ^ = ZZ ,( - 8 ^ as before (with /' = 0). 
However, since e = 1, we may also take f' = l, and we obtain ZZ^ 9 - 1 + JL^ = ZZ^ 8 - 1 * + ZZ® = 
(x 5 + x 4 + x 2 + 1, x 3 + x + 1) G MR'(s) too. This is equivalent to taking /' = x + 1 in 
Corollary [9T4l 

Remarks 9.10 (i) We could replace // 6y any element o/MP(s'). (m) i?y construction, 
H* — g/i — A • /x' where q = A' • x e + /' and deg(g) = e. 
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9.3 A Characterisation 



In this subsection, we give an analogue of Theorem 15.91 for Ann*(s) when /i is unique 
and fiQ = 0: we characterise {/ G Ann*(s) : deg(/) is minimal} in this case. This 
was stated without proof for D = ¥ in [181 proof of Theorem 3.7(iii)]. We also describe 
h( x ) = S^= S q i (/ ' §)j xi '■ (The case e < is an easy consequence of Theorem 15.91 ) 
First a simple but useful proposition. 

Proposition 9.11 Let f,r E D[x) and n — 1 > k = deg(f) — deg(r) > 1. If f, f + r E 

Ann(s) then r E Ann(s^ n ~ fc ^). 

Proof. As n — k > 1, Ann(s^ n_fc ^) is well-defined. Let d = deg(/ + r) = deg(/) = 
deg(r) + k and j satisfy deg(r) — n + k < j < — 1. Also, d — n = deg(r) — (n — k) and so 
(r • s)j = ((/ + r) • s)j — (/ ■ s)j = and hence r E Ann^ 71- ^) . ■ 

It is convenient to have the notion of a 'jump point'. 

Definition 9.12 (Jump point) Let n > 2, s E D n and 2 < j < n. We say that j is a 
jump point of s ifLj > Lj_i and write J(s) for the set of jump points of s. 

We do not assume that J(s) ^ 0. Evidently, the following are equivalent: (i) j E J(s) (ii) 
Aj 7^ and e^-i > (iii) j' = j — 1 (iv) Lj = j — Lj_i > Lj. In Theorem 19.91 s n+ \ was 
chosen so that n + 1 is a jump point of (si, . . . , s n+ i). 

Theorem 9.13 (Cf. fJM Theorem 3.7(m)]) Let n > 2, s E D n , fi = 0, e > and 
£ = /iL +1 - Then f E Ann*(s) has minimal degree if and only if 

deg(/) = n + 1 — L, n' > 1 and £ • f = q-Jt + f 

where q — (£ ■ f) div /j,, r — (£ ■ /) mod G MP(s') and r 7^ 0. In particular, if D is a 
field then f = q ■ fi + r where deg(g) = e and r — a - ~p' for some a E D x . 

Proof. Corollary 19.41 implies that qjl + r E MR'(s) and hence so is /. Conversely, let 
/ be as stated. Then deg(/) = n + 1 — L ra by Corollary 19.51 Since /i = 0, L > 1 and 
e > implies that L < |. Pseudo-dividing £ ■ f by \i (as in [3 Algorithm R, p. 407]) 
gives £ ■ f = qji + r with deg(g) = n + 1 — 2L = e>0 and r = or deg(r) < L. 
Since r = £fo 7^ 0, < d = deg(r) < L. Let k — L — d > 1. Then n — k > 1 since 
n>2L>L>L — d = k and s( n_fe ) is well-defined. By Proposition 19. Ill r G Ann(s^ n ~ fc ^). 

First we show that n — k = n'. Since d < L and n' is maximal by [121 Proposition 4.1], 
n — k < n' . In particular, n' > 1 and we are done if n' = 1. Suppose that n — k < n' . If 
L' < d+1 then 

n-L + L'<n-L + d+ l = n- A; + l<n / = L + L / -l (2) 

and e < 0. Similarly, L' = d is impossible and so d < L'. To treat this case, we generalize 
fl2]) by partitioning [Li,L') into [L 1; L^ ),..., [L Jp l , L Jp ) for some integer p. There is a 
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unique integer t > 2 with d £ [Lj t _ x , Lj t ), where j tl = 1. Since d < Lj t , we must have 
n — k < j t — I a.s before and therefore 

n — L + hj t _ 1 < n — L + d = n — k < j t — 1 = Lj t + L J - t _ 1 — 1 < L + L J - t _ 1 — 1 

which implies that e < 0. We conclude that n — k = n' and L' < d as r £ Ann(s'). If 
L' < d, then (n-L) + L' + l<(n-L) + d = n- £; = n / = L + L / -l and e < -1, so 
that d = L'. 

Finally, if D is a field, we have shown that / = q ■ fi + r where r £ MP(s') and ro ^ 0. 
Moreover, L' < L — 1 = n' — L' so that e n > = n' + 1 — 2L' > 1 and // is the unique monic 
minimal polynomial of s' by [T2l Corollary 3.27 or Theorem 4.16]. Hence r = a fx' for some 
aeD x . ■ 



Remark 9.14 In Theorem \9.13[ we have deg(g) + deg(/i) = ra + l — L<n since L > 1 
and follows easily from fl^ Corollary 3.25] (which does not require D to be a field), 
that £ ■ f = (g/i + r, (q/j, + r) 2 ) = qjt + r. 

Example !5.11l (cont.) From Theorem 19. 131 and Remark 19.141 

MR'(s) = {(x 5 + x + 1, x 3 + x 2 ), (x 5 + x 4 + x 2 + 1, x 3 + x + 1)}. 

Appendix 

We give an alternative proof of Lemma 19.21 using reciprocals of annihilators and apply it 
to the complexity of reversed sequences. This generalises Theorem 3 of jU p. 149] (which 
was proved using using Hankel matrices over a field) to a factorial domain. 

Let s denote the reverse of s £ D n and L = L(s). We assume that (j,^ £ MP(s). We 
begin with a basic result on annihilators and reciprocals. 

Proposition 9.15 Let f £ Ann(s) and d = deg(f) and let k > be the highest power of 
x dividing f . Then f* £ Ann(s n , . . . , sjt+i). In particular, if x \ f then f* £ Ann(s). 

Proof. First assume that k = 0. For d — n < j < — 1, 

= (f(x- l )s)^ d = (/ • ( Sn X + • • • + Sl X n )) d -j = (X n+1 fs) d ^ = (J^^-n-i = 

since d — n < d — j — n — 1 < — 1. In general, let /' = Then /' £ Ann(sfc+i, . . . , s n ) 

by Proposition I9.1( i) and /* = /'* £ Ann(s n , . . . , Sk+i) by the case k = 0. ■ 

Next we give an alternative proof of Lemma 19.21 using reciprocals of annihilators. 

Lemma 9.16 Let n > 2, f,g £ Ann(s). If x\f, f/x £" Ann(s^ n_1 - ) ) and x \ g, then 
deg{g) >n + l-deg(/). 
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Proof. Since go 7^ 0, g ^ and g* G Ann(s) by Proposition 19. 15( ii) . Let k be the highest 
power of x dividing / G Ann(s) and / = f/x k , so that deg(/*) = deg(/) = deg(/) — k. 
Firstly, / G Ann(sfc + i, . . . , s n ) by Proposition I9.15( i) and x\ f, so /* G Ann(s n , . . . , s^+i) 
by Proposition I9.15( ii). But /* ^ Ann(s n , . . . , S&), for otherwise / G Ann(sfc, . . . , s n ) and 
hence //a; = x k ~ 1 f G Ann(s) by Proposition 19.151 This is impossible since Ann(s) C 
Ann(s ( ^ n_1 - ) ). We can now apply Lemma [3.31 to /* and g*, yielding deg(g*) > n — k + 1 — 
deg(/*) = n + 1 - deg(/) i.e. deg(^) > n + 1 - deg(/). ■ 



Lemma 9.17 (i) If fio ^ then L < L and if fio = then L > n + 1 — L. 

Proof. If fi ^ then /i* G Ann(s) by Proposition 19. 15( ii) . so L < deg(/i*) = deg(/z) = L. 
Suppose now that (1q = and let x fc |/i with > 1 maximal and put g = [ijx k . Then 
go ^ 0, deg(g) = L — and (7 G Ann(sfc + i, . . . , s n ) by Proposition 19.15( h). Hence g* G 
Ann(s n , . . . , Sk+i). However, g* £ Ann(sfc, . . . , s n ), for otherwise g G Ann(sfc, . . . , s n ) and 
fi/x = x k ~ l g G Ann(s) which contradicts the minimality of \x. We now apply Lemma [3.31 

to 9*- 

L > L n _ fc _i >n — k + l — deg(g*) = n + 1 - L. 



Example 9.18 Lei D = F 2 and s = (1, 1,0,0). One easily checks that x 2 is the unique 
minimal polynomial. By Lemma 9.17\ L > 3. In fact // 3 ) = x 3 + x 2 and so L = 3. Thus 



the bound of Lemma 9.17 is tight 



We will now show that in the special case n = 2L that L = L + 1 when D is factorial 
and jj Q n ^ = 0. Recall that if D is a factorial domain and n > 2L then s has a unique 
minimal polynomial. 

Corollary 9.19 (Cf. f^j Theorem 3]). Let D be a factorial domain, n = 2L and fi be the 

unique minimal polynomial of s . Then 



L 



L i/ fi ^ 
L + 1 if fi = 0. 



Proof, (i) We have L < L by Lemma ET7(i). Let /2 G MP(S). If /a ^ 0, /2* G Ann(s), 

L = L < L and we are done. If /io = 0, Lemma [9.17( h) implies that L > n + 1 — L = 
2L + 1 - L. Hence L > L + 1 > L and L = L. 

(ii) Suppose that /io = 0. Let / be such that / G Ann(s), /o / and deg(/) is 
minimal. Then deg(/) = n + 1 — L by Lemma 19.21 and L < n + 1 — L<L + lby 
Proposition 19.151 and the fact that n = 2L. We now show that L < L is impossible. 

Let / e MP(S). If f ^ 0, then /* G Ann(s) and so^L < deg(/) = L. We cannot 
have L = L, for then f* = fi and /io 7^ 0. Hence L < L and we cannot have L < L. 
By the first paragraph, L = L + 1. On the other hand, if fo = 0, then by Lemma [9.171 
L>n + 1 — L = L + 1. Thus in either case, L > L + 1 and so L = L + 1. ■ 
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Typographical Error 

In [10], immediately after Definition 3.3, A(/, s) = Ylt=o fk s n-d+k- 
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